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ABSTRACT 

Thanks to the Dirac equation, the hydrogen-like atom at high Z offers a 
precise model of relativistic bound state, allowing to test properties of un¬ 
polarized and polarized structure functions analogous to the hadronic ones, 
in particular: Sivers effect, sum rules for the vector, axial, tensor charges 
and for the magnetic moment, positivity constraints, sea contributions and 
fracture functions. 
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1 INTRODUCTION 

In this work we will study the hydrogen-like atom, treated by the Dirac 
equation, as a precise model of relativistic two-particle bound states when 
one of the constituent is very massive. We consider the case of large Z so 
that Za ~ 1 and relativistic effects are enhanced. What we neglect is 

• the nuclear recoil 

• the nuclear spin (but not necessarily the nuclear size) 

• radiatives corrections, e.g. the Lamb shift ~ a{ZaYme- 
In analogy with the quark distributions, we will study: 

• the electron densities q(kT,k~^) or q(h,k^) where b = {x,y) is 

the impact parameter; 

• the corresponding polarized densities; 

• the sum rules for the vector, axial and tensor charges and for the atom 
magnetic moment; 

• the positivity constraints; 

• the electron - positron sea and the fracture functions. 

We use k^ = k^ + kz instead of the Bjorken scaling variable k^/Patom which 
would be of order nie/Matom, hence very small. \k~^\ can run up to Matomi 
but typically \k'^ — m\ ~ Zam. 

2 Deep-inelastic probes of the electron state 

Deep-inelastic reactions on the atom are, for instance: 

• Compton scattering: y* -f- e“(bound) —7/ + ej(free), 

• Moeller or Bhabha scattering (replacing the 7 by a e^), 

• annihilation: e"*"-|-e“ (bound) —7-|-7, 
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the Mandelstam invariants s, t and u being large compared to rrig. In Comp¬ 
ton scattering for instance, taking the z axis along Q = k^ —k/, the particles 
7i, 7/ and ej move almost in the —2; direction and we have 

( 1 ) 

k-r ~ kjj’ + kj^ — kjj. = —PfT{nudeus). (2) 


3 Joint (b, k^) distribution 


In the ’’inhnite momentnm frame” Pz S> M, deep inelastic scattering mea- 
snres the gange-invariant mechanical longitndinal momentnm 


kz — {Pz/Matom) frame ~ Pz -^z ( 3 ) 

and not the canonical one Pz = —id/dz- It is not possible to dehne a joint 
distribntion q{)^Ti kz) in a gange invariant way, since the three qnantnm op¬ 
erators ki = —idi — Ai{t, X, y, z) are not all commnting. On the other hand we 
can dehne nnambignonsly the joint distribntion g(b, kz). Note that q{h, 
can be measnred in atom -|- atom collisions where two hard snb-collisions 
occnr simnltaneonsly: {ej” -|- 62 and Ni -|- iV2} or: (ej” -|- N2 and -|- A^i}. 
Given the Dirac wave fnnction \&(r) at t = 0 in the atom frame, we have 


dNe-/[d^h dk+/{ 27 r)] = q{h,k+) = fc+) $(b, fc+) , 

<h(b, A;”'”) = J dz exp - 1 -— ix(b, 2:)| <h(r), 


$(r) = 


d'i(r) -h d' 3 (r) 
d'ofrj - tkdfrj 


x{h, z) = / dz' V{x,y,z') = —Za 

J ZQ 


sinh 


-1 


— sinh ^ 


-1 /^_o 
b 


( 4 ) 

( 5 ) 

( 6 ) 
( 7 ) 


The two-component spinor $ represents The ’’gange link” exp{—ix(b, z)} 

transforms T in the Conlomb gange to T in the gange A'^ = 0 {Az = 0 in 
the inhnite momentnm frame). The choice of zq corresponds to a residnal 
gange freedom. 

One can check that Q is invariant nnder the gange transformation V (r) —> 
l/(r)-|- Constant, Pn En+ Constant. Snch a transformation may resnlt 
from the addition of electrons in onter shells; it does not change the mechan¬ 
ical 4 -momentnm of a K-shell electron. 
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4 Joint (ky, k^) distribution 

The amplitude of the Compton reaction is given, modulo a matrices, by 

oc J dh e-*Q--*k/ r-*x(b,.) _ $(kr,A;+) (8) 

with zq = —oo. exp{—ikj ■ r — iy} is the final wave function distorted by 
the Coulomb potential, in the eikonal approximation. Equivalently, 

$(k-r, A;+) = J dh ^-ikr-h-ik+z-ixir) _ J ^2^^ g-ik^-b ^ 

with the identification = En + kf^z + Qz- [check: in a semi-classical 
approach, fc°(r) = En — C(r) and kf^z + Qz = kz{r) — E(r) at the collision 
point]. For the annihilation reaction, the incoming positron wave is distorted, 
then zq = -|-oo. Thus, the gauge link takes into account either an initial or a 
final state interaction mum- The quantity 

q{kT, k+) = <f)t(kT, k+) <f)(kT, k+) (10) 

depends on the deep inelastic probe, contrary to g(b, k~^) and 

g(fc+) = J g(b, k^) d% = J q(kT, k^) d^k^/ (27rf . (11) 

5 Spin dependence of the electron density 

Ignoring nuclear spin, the atom spin is j = L + s. We will consider a j = 
1/2 state and denote by = 2(j) and S® = 2(s) the atom and electron 
polarization vectors. and S® without bar indicate pure spin states. The 
unpolarized electron density in (b, k~^) space from a polarized atom is 

g(b,A;+;S^) = <l>^(b,A;+;S^) <l)(b, A;+; S^) (12) 

and the electron polarisation is given by 

S"(b,fc+;S^) q{h,k+-,S^) = $^(b,A:+;S^) a $(b,fc+;S^). (13) 

Taking into account the conservations of parity and angular momentum, 
the fully polarized density can be written as 

g(b, fc+, S^) = qib, fc+) [1 + Con (S^ ■ n) + (S^ ■ n) + Cnn (S^ ■ n) 
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+ Cu StS^ + Ci^ StiS^ ■ n) + C^i (S^ • if) (S^ • ^)(S^ ■ tt)] (14) 

where -ff = h/b and h = z x d. The Qj-’s are functions of b and k~^. 

Similar equations work for instead of b. The link with Amsterdam 


notations |3] is, omitting kinematical factors, 

q{kT, k^) = fi fi Con = fiT (15) 

/i Cu = gi fi Cno = -h^ (16) 

fi Cnn = hi- fi Ci-^ = giT (17) 

/i C'^TT = hi + hiT fi C^i = hiL , (18) 

The b - or kr integration washes out all correlations but Cu and Ctt = 
1 [Cnn + Ctttt], giving 

g(fc+,S";S^) = g(A;+) + Aq{k^) + 6q{k+) . (19) 


5.1 Sum rules 

Integrating P|) over one obtains the vector, axial and tensor charges 

q= j q{k'^) dk'^ / {2n) = / T^(r; S"^) T(r; any S^) (20) 

Ag = y Ag(A;+) dA;+/(27r) = J dh ^^(r; S^) T(r; = z) (21) 

Sq = J 6q{k+) dk+ /{ 27 r) = J d^r ^^(r; S^) fd T(r; = x). (22) 

For the hydrogen ground state, 

g = l, Ag=(l-eV3)/(l + a, dg = (1 + eV3)/(l + ^ (23) 

where ^ = Za/{1 + 7), 7 = E/me = y^l — (Za)^. 

Note the large ’’spin crisis” Ag = 1/3 for Za = 1. 

5.2 Results for the polarized densities in (b,/;:^) 

The 2-component wave functions of a A = -1-1/2 state write 
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Other orientations of can be obtained by rotation in spinor space. 

For the (b, k'^) representation we can ignore the second term of ((7j). Then 
n(b, k'^) and w(b, k^) are real and given by 


V 

w 



& \ 

r + i^z ) 


^-ik+z+iEz-ix{h,z) f(^f)/r 


(25) 


where /(r) cx r'’' ^ exp(—mZar) is the radial wave function. Then, 


q{b, k+) 
Cnn{b,k^) 
Conib, k^)=Cnoib, 

Caib,k+) = C^^ib,k+) 
CUb,k^)=C^iib,k^) 


W + V 
1 

—2 wu /+ v^) 

+ v'^) 

0 . 


(26) 


5.2.1 Sum rule for the atom magnetic moment 

Consider a classical object at rest, of mass M, charge Q, spin J and time- 
averaged magnetic moment fl. Its center of mass r^ and the average center 
of charge (re) coincide, say at r = 0. Upon a boost of velocity v, the center 
of energy r^ and (rc) are displaced laterally by 

bG = vxJ/M, (be) = V X/I/Q . (27) 

be and (be) coincide if the gyromagnetic ratio has the Dirac value Q/M. For 
the hydrogen atom, bg is negligible and the magnetic moment is almost fully 
anomalous. In the inhnite momentum frame (v ~ z), we have an electric 
dipole moment [H] 

— e (b) = p-A z X , (28) 

the transverse asymetry of the b distribution coming from the Con term of 
(HI. We recover the atom magnetic moment 

=-e (1 + 27)/(6me) (29) 

(the anomalous magnetic moment of the electron being omitted). 
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5.3 Results for the polarized densities in (ky, k~^) 

For the (k^, representation we should take zo = ±cx) but then ((71) di¬ 
verges. In practice we assume some screening of the Coulomb potential and 
take \zo\ large but hnite, which gives 


X(b,^) 


—Za 


±ln(2|2:o|/&) + sinh ^{z/h) , 


(30) 


the upper sign corresponding to Compton scattering and the lower sign to 
annihilation. Modulo an overall phase, 



POO 

2tt bdb 6^*^" 
Jo 


/ Joikrb) w{b, k+) \ 
\MkTb)v{b,k+)) ’ 


(31) 


q{kT, k+) 
Cnni^kj ', k ) 
Con{kT,k+) = Cr,o{kT,k+) 

Cu{kT,k'^)=CUkT,k+) 
G.(fcT,fc+) = -aKfcT,fc+) 


I ~ |2 I I ~|2 

\w\ + \v\ 

1 

2Q^(h*ti})/(|u}p -f- |hp) 

(Itcp - \v\^)/{\w\^ -f |h|^) 

2gfJ(rh})/(|M;p-F IhH. (32) 


The factor 5“*^" (Compton case) behaves like a converging cylindrical 
wave. Multiplying $(r), it mimics an additional momentum of the electron 
toward the 2; axis. In fact it takes into account the ’’focusing” of the hnal 
particle by the Coulomb held [S]. It also provides the relative phase between 
w and V which gives non-zero C'on(fcr, k~^) and Cno{kT, k~^) (Sivers and Boer- 
Mulders-Tangerman effects). Similarly, (annihilation case) takes into 

account the defocusing of the incoming positron. 


5.4 Positivity constraints 

The spin correlations between the electron and the atom can be encoded in a 
” grand density matrix” R [H] , which is the hnal density matrix of the crossed 
reaction nucleus atom{S^) + e’''(—S®). 

Besides the trivial conditions \Cij\ < 1 the positivity of R gives 

(1 ± Cr^nf > [CnO ± C,r.f + {Ql ± + [C^l T (33) 

These two inequalities as well as |Ci;| < 1 are saturated by (l7F)|l and (IH^ . In 
fact R is found to be of rank one. It means that the spin information of the 
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atom is fully transferred to the electron, once the other degrees of freedom 
and b or k-p) have been hxed. If there is additional electrons or if we 
integrate over for instance, some information is lost and some positivity 
conditions get non-saturated. Conversely, the hypothesis that R is of rank 
one leaves only two possibilities: 

Cnn = ±1 , Con = ±^^0 , , (34) 

with saturated. The hydrogen ground state choses the upper sign. After 
integration over b or kj-, we are left with the Soffer inequality, 

2 |5g(fc+)| < g(fc+) + Ag(A:+) (35) 

which in our case is saturated, even after integration over see (j2Sl)- 

6 The electron-positron sea 


The charge rule (EOl) involves positive contributions of both positive and 
negative values of . So it seems that there is less than one electron (with 
k'^ > 0) in the atom. This paradox is solved by the introduction of the 
electron-positron sea. 

Let \n) be an electron state in the Coulomb held. Negative n’s are as¬ 
signed to negative energy scattering states. Positive n’s up to ns label the 
bound states (—m < < +m) and the remaining ones from -|- 1 to 

-|-cx) are assigned to unbound states of positive energy, En > +m, considered 
as discrete. Let \k,s) be the plane wave of four-momentum k and spin s, 
solution the free Dirac equation. The destruction and creation operators in 
these two bases are related by 

ak,s = ^ {k,s\n) an, = ^ ai{n\k,s). (36) 

n n 

In the Dirac hole theory, the hydrogen-like atom is in the Fock state 

\Hn) = ttn all (^-2 ■ ■ ■ ®loo |QED-bare nucleus) (37) 

The number of electrons in the state \k, s) is 

A^'^'^(k,s) = (4,, ak,s) = \ {k,s\n)\^+ ^ \{k,s\n')\^ 

n'<0 
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(38) 


For a nucleus alone (but ” QED-dressed”) the first factor of (pTTfl is missing 
and the hrst term of is absent. Therefore, passing to the continuum limit 
{k,s\n) —> $(kT,/c+), 

r 

9ir 

One term n' of (jd8j) corresponds, e.g., to Compton scattering on an elec¬ 
tron of the Dirac sea, producing a fast electron plus the positron \n') = 
vacant. It is the fracture function of reaction 'ji + A^'jf + A + e'^ + X. 
Similarly the number of positrons is (with > 0): 

]V>”“‘(k.s) = {o_,,_,al,,_,>= i: |{-4,-sK)|^ (40) 

0<n'^n 

For a nucleus alone the condition n' ^ n is relaxed. By difference, 

NSSL - 'vrr = i ^ 9(V) (41) 

Jk+<o In 

One term of Eq. iOD corresponds to the extraction of an electron of large 
negative energy, giving a fast positron and an electron in \n'). Eqs. fP|) + 
m and (EOD tell that the atom and nuclear charges differ by e. But the 
charge of the electron-positron virtual cloud surrounding the nucleus is not 
zero (charge renormalisation). 

7 CONCLUSION 

We have seen that the leading twist structure functions of the hydrogen¬ 
like atom at large Z has many properties that are supposed or verihed for 
the hadronic ones, in particular: sum rules, longitudinal ’’spin crisis”, Sivers 
effect, transverse electric dipole moment in the Poo frame, etc. It remains 
to evaluate these effects quantitatively. With this ’’theoretical laboratory” 
one may also investigate spin effects in fracture functions, non-leading twist 
structure functions, Isgur-Wise form factors, etc. The electron-positron sea 
may deserve further studies: to what extent is it polarized or asymmetrical 
in charge ? Is the charge renormalisation of the nucleus found in Sec.6 the 
same as in standard QED ? Is it hnite and calculable for an extended nucleus 
? 

Finally it may be interesting to do or re-analyse experiments on deep 
inelastic Compton scattering (’’Compton prohle” measurements). 
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